Instability of the origami of a ferrofluid drop in a magnetic field 
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Capillary origami is the wrapping of an usual fluid drop by a planar elastic membrane due to the 
interplay between capillary and elastic forces. Here, we use a drop of magnetic fluid whose shape 
is known to strongly depend on an applied magnetic field. We study the quasi-static and dynam- 
ical behaviors of such a magnetic capillary origami. We report the observation of an overturning 
instability that the origami undergoes at a critical magnetic field. This instability is triggered by 
an interplay between magnetic and gravitational energies in agreement with the theory presented 
here. Additional effects of elasticity and capillarity on this instability are also discussed. 

PACS numbers: 47.65.Cb,46.32.-|-x,68.08.-p 



Generally, a solid in contact with a static liquid in- 
terface is undeformed by surface tension forces at large 
scales. However, for sub-millimetric scales or when grav- 
ity is negligible, capillary forces may deform an elas- 
tic structure in various domains (see [H for a review): 
adhesion (coalesence of wet hairs Q), biological sys- 
tems (floating flowers 0), or industrial applications (mi- 
croscale fabrication such as 3D photovoltaic cells or 
MEMS [H ) . The most well-known phenomenon involving 
the interplay between elasticity and capillarity is capil- 
lary origami. It concerns the spontaneous wrapping of a 
liquid droplet by a planar elastic sheet when the capillary 
forces dominate the restoring elastic ones 0, 0] . Various 
folding shapes (spherical, cubic or triangular encapsula- 
tion) can then be tuned from the geometry of the initial 
flat membrane. A new challenge is to accurately control 
the folding and unfolding of the origami without the use 
of mechanical part. Such a capillary origami control has 
been recently performed by means of an electric field and 
could lead to potential applications to digital displays 
@,01- However, the capillary origami is generally consid- 
ered as quasi-static phenomenon and dynamical studies 
are rare. To our knowledge, only one experiment con- 
cerns the dynamical selection of the final shape of the 
capillary origami when a drop impacts an elastic mem- 
brane [8|. 

In this Letter, a drop of magnetic fiuid is deposited 
on a flat elastic membrane and is submitted to a mag- 
netic field. The shape of the ferrofiuid drop is known 
to strongly depend on the applied magnetic field 0-11 



The wrapping of the drop by the thin elastic membrane 
(origami) is thus expected to be strongly modified by the 
magnetic field. Both quasi-static and dynamical behav- 
iors of such a magnetic capillary origami are reported 
here. The most striking one is the observation of an 
overturning instability that the origami undergoes at a 
critical magnetic field. This instability is shown to arise 
from the interplay between magnetic and gravitational 
energies. 

The ferrofluid used is an ionic aqueous suspension syn- 
thesized with 12.4% by volume of maghemite particles 



(Fe203 ; 7 ±0.3 nm in diameter) [l2|. The properties of 
this magnetic fluid are: density p = 1550 kg/m'^, surface 
tension 7 = 43 ± 3 mN/m, initial magnetic susceptibility 
Xi = 0.75, magnetic saturation M^at = 36 x 10'^ A/m, 
and dynamic viscosity 1.4 xlO~'^ N s/m^. The elastic 
membranes arc made of polydimethylsiloxane (PDMS - 
Dow Corning Sylgard 184), a 10:1 polymer/curing agent 
mix. The PDMS is spin coated at rotation rates in the 
range 1800 - 2800 rpm over Emery polishing paper (grit 
0) with average roughness of 192 fim. The thickness h of 
the membrane depends on the rotation rate and ranges 
from 50 to 100 /im. Its roughness significantly reduces 
the adhesion of the membrane on the substrate Q. The 
ferrofluid drop alone or wrapped by the elastic membrane 
(origami) is then placed between two horizontal coils, 25 
cm (rcsp. 53 cm) in inner (resp. in outer) diameter, 7 cm 
far apart. A dc current, /, is supplied to the coils in series 
by a power supply (50 V/50 A). The vertical magnetic 
induction, B, generated is up to 780 G and is measured 
by a Hall probe located near the drop. The magnitude of 
B is proportional to / and is controlled cither manually 
or by means of a ramp generator of typical speed rate 
0.4 A/s. Deformations of the drop are visualized from 
a side view by means of a 45 degree mirror reflecting a 
diffuse lighting towards a high-resolution camera (Pix- 
elink 2208 x 3000 pixels) located above the drop, and are 
recorded up to a 150 Hz sampling. 

We first study the deformation of a ferrofluid drop 
alone deposited on a superhydrophobic surface and sub- 
mitted to a vertical magnetic field. As shown in the inset 
of Fig. [I] the drop is quasi-spherical at B = (a). When 
B is increased, the drop lengthens in the direction of B 
and its shape changes continuously from a semi-cUipsoid 
(b), then a pointed ellipsoid (c) up to a sharp tip (d). 
By analogy with a semi-ellipsoid, we define a (resp. h) 
as the semi-major (resp. minor) axis of the deformed 
drop, a is measured from the top of the drop to the 
horizontal plane where its width is largest, whereas h is 
its half-maximal width (see inset of Fig. [1]). Both are 
measured on the image of the drop with a 5% accuracy. 
The aspect ratio a/h of the drop is plotted in Fig. [l]as 
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FIG. 1. (Color online) Aspect ratio a/fe of a ferrofluid drop 
on a superhydrophobic substrate versus the magnetic Bond 
number Bom, for two values of the initial drop radius: _Ro = 

I. 2 (A) and 0.93 mm (□). Inset: Photo of the ferrofluid drop 
deformation for different vertical magnetic fields: B = (a), 
260 (b), 520 (c), and 780 G (d) corresponding to Bo^ = 0, 

II, 31, and 53 respectively. -Ro = 1.2 mm. Bar scale is 2 mm. 



a function of the magnetic Bond number, Boi„. This di- 
mensionless number characterizes the order of magnitude 
of the ratio between magnetic and capillary energies . 
Bo„i = x{B)B'^ Rq / {p.Q"f) with i?o the initial drop radius 
at i? = 0, /iQ = 47r X 10^' H/m the magnetic permeabil- 
ity of the vacuum, and x(-B) the magnetic susceptibility 
of the ferrofluid which is a known decreasing function of 
B with Xi = x{B = 0) [l^l. The elongation of the drop is 
found to strongly increase with Bom since the ferrofluid 
tends to align towards the direction of B. Moreover, as 
shown in Fig. [1] both curves performed for two different 
i?o superimpose underlying that a small drop is less de- 
formed than a larger one for the same applied B. Note 
that such a deformation has been also observed when the 
ferrofluid drop is surrounded by a non misciblc fluid of al- 
most same density leading to a full ellipsoidal shape: a/b 
is then found to increase continuously with B for small 
Xi 1 whereas for large one (xi ^ 20) a discontinuous 
deformation occurs [l^l. The elongation of a full ellip- 
soidal droplet can be described theoretically in the limit 
of a linear [13, [13 or nonlinear drop magnetization. 
To our knowledge, no analogous analytical computation 
exists for a ferrofluid drop with a semi-ellipsoid shape (as 
the one in Fig. [T]). 

Let us now focus on the behavior of a ferrofluid drop 
wrapped by an elastic membrane and submitted to a 
magnetic field. A geometric shape of the membrane is 
manually cut out from the PDMS layer and placed on 
the superhydrophobic surface. The membrane shape is 
chosen to be an equilateral triangle of side L (5 < L < 15 
mm) such that the closed state of the origami will have 
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FIG. 2. Deformation and instability of the magnetic capillary 
origami as the vertical magnetic field is increased: _B = (a), 
150 (b), 400 (c), and 600 G (d) corresponding to Bo,„ = 0, 
10, 48, and 84 respectively. An overturning instability occurs 
between pictures (c) and (d) (see text and movie in Ref. [isl]'). 
i?o = 2.7 mm. L = 10 mm. /i = 65 /im. Bar scale is 2 mm. 



a pyramid- like shape [J]. A typical experiment is as fol- 
lows. First, at i? = 0, a ferrofluid drop is deposited on 
the planar membrane. Second, the volume of the drop is 
adjusted such that, due to the competition between capil- 
lary and elastic forces, the three corners of the membrane 
touch themselves and wrap the spherical drop as shown in 
Fig. When B is turned on and is increased, the drop 
lengthens vertically (Fig. [2Jd), then exhibits a pyramidal 
shape at higher B due to the presence of the membrane 
(see Fig. For a critical value B^, corresponding to a 
critical Bond, BoJ^,-,, the origami suddenly overturns, then 
oscillates around its vertical position before stopping its 
rocking motion (see movie in [15[). This leads to a new 
static configuration of the origami (see Fig.[2}i): the drop 
now forms a cone at its top, and rests on one corner of 
the membrane, the two other corners wrapping it hori- 
zontally. If B is now decreased, the origami falls forward, 
and never returns to its initial configuration thus showing 
a hysteretic behavior. 

To quantitatively study the above instability, the max- 
imum height Zmax of the magnetic capillary origami is 
measured during a linear temporal increase of _B at a 
tunable rate ranging from 3 G/s to 40 G/s. This rate 
is slow enough to consider a quasi-static evolution of B 
(see below). Figure [3] shows the dimensionless height 
ZmaxI Ro as a function of Bom for the origami and for 
a ferrofluid drop alone (i?o being the initial radius). For 
the drop alone, z^ax/Ro is a continuous function of BOm. 
For the origami, a jump of z^ax/Ro is observed for a crit- 
ical Bond number, Bo^j^, corresponding to the threshold 
of the overturning instability for which the origami shape 
changes from a pyramid (Fig. [2j:) to a cone-like shape 
(Fig. [2li). Due to inertia, this jump in Zmax/Ro is accom- 
panied by an overshoot {zmax/Ro reaching a maximum 
higher than the equilibrium) followed by underdarnped 
oscillations (see bottom inset of Fig.[3]and movie in jl5|). 
Moreover, just after the instability occurring at BoJ^^, the 
origami in its new configuration has a height close to the 
one of the drop alone meaning that the stress applied by 
the membrane has relaxed. When iterating this exper- 
iment for different sizes L of the triangular membrane 
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FIG. 3. (Color online) Height of the origami [red (light gray) 
dots] rescaled by its initial radius {Ro = 2.3 mm), as a func- 
tion of the magnetic Bond number Bom, and the correspond- 
ing evolution for a ferrofluid drop alone (black dots, Rq = 2 
mm). L — 8.4 mm. h — 56 fim. A jump in the height of 
the origami is observed for a critical Bond BoJii — 25.3, cor- 
responding to the instability shown between Fig. [2}; and 2d. 
Bottom inset: zoom near the instability region. Top inset: 
Bon, versus the size L of the elastic membrane. 



of fixed thickness h, the onset of this instability Bojjj is 
found to increase with L, as shown in the top inset of 
Fig.E 

We first explain these above observations by dimen- 
sional scaling laws. At the onset of the instability, an 
increase of the origami height is observed (see Fig. [3]) un- 
derlying that gravitational energy Eg should be taken 
into account. One has Eg ~ pVgZmc with z,„c the 
mass center of the origami and V the volume of the 
ferrofluid. One can assume ^ L, and V ^ 
since the initial configuration of the origami is quasi- 
spherical, and thus Eg ~ L*. The magnetic energy writes 
Em ~ pqMHV ^ L^B'^x{B)/ Ho since the magnetiza- 
tion M ~ x{H)H and the magnetic field H = B/^q 
[Tst . The balance between the magnetic and gravita- 
tional energies thus leads to the critical magnetic field 
Be at the instability threshold such that x{Bc)B'^ ^ L. 
Thus, Bojjj ^ x{Bc)B1 is expected to increase with L 
as observed in the inset of Fig. |31 Bojj^ is also expected 
to be independent of the membrane thickness h as found 
experimentally. A similar scaling law analysis balancing 
the magnetic energy with the capillary energy, Ec, or the 
elastic one, E^, leads for both cases to inverse predic- 
tions, i.e. a decrease of BoJ^,-, with L since Ec ~ 7-^^, 
and Ee ^ Eh^K^L"^ {E being the Young modulus of the 
membrane of curvature k ^/L). Thus, the origami 
overturning occurring at Be is triggered by a competition 
between magnetic and gravitational energies. Note that 
the role of capillarity and elasticity is non zero. When 
-B = 0, it drives the origami formation when Ec > E^ 
(i.e. for L > ^ Eh? [4|), and they both deform the 
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FIG. 4. (Color online) Variation of the height of the origami 
mass center, Az^c, during the instability as a function of 
/(, AMr„{H)dHc for different sizes L of the membrane: L = 
7.8 (•), 8.5 (■), 9.4 (t), and 10 mm (♦). Membrane thickness 
ft = 65 ^m. Solid line has a slope ^o/{2pg). 



fcrrofiuid shape when B ^ 0. 

Let us now characterize more quantitatively this in- 
stability. To do that, one has to compute the mag- 
netic energy E^ = -^^oV M{H)dHo with M the fer- 
rofluid magnetization assumed uniform within the drop 
14| . H is the magnetic field within the ferrofluid droplet 
that depends on the applied magnetic field Hq through 
the implicit equation H = Ho — DAI{H), D being 
the demagnetization coefficient that depends on the fer- 
rofluid shape For a full ellipsoid, it is known that 
[In — 2e] with e = \J\ — (b/a)"^ the ecccn- 
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tricity, and a (resp. b) the semi-major (resp. minor) axis 
lol 11 [. For instance, in the case of an ellipsoidal drop 
alone, = for a drop infinitely stretched in the direc- 
tion of Hq, Z? = 1/3 for a spherical drop, and D = 1 for 
an infinitely flat drop normal to Hq . The magnetic energy 
Em thus depends on the external magnetic field Hq and 
on the droplet shape through its magnetization M (H) 
since iZ is a function of both D and Hq. Consequently, 
the computation of Em demands to integrate M{H) over 
the external magnetic field Hq. We estimate A4{H) 

by using the usual Langevin's expression: M{H) = 

Msat^3x^H/Msat) whcrc Cix) = coth(a^) - 1/x m. 
The initial magnetic susceptibility, Xi = ^-\h^o, and 
the magnetic saturation, Msat = M{H)\h^oo, are given 
by the ferrofiuid properties (see above). 

The variation of the magnetic energy during the in- 
stability, occurring at the critical field Hq = He, reads 
AEm = -^^aV j"' AM{H)dHe, where AM = M{H2) - 
M{Hi) is the variation of M before and after the in- 
stability (denoted by subscripts 1 and 2 respectively). 
Hi= He- DiM{Hi) and H2 ^ H^ - D2M{H2) are the 
magnetic field within the drop before and after the insta- 
bility. Since the origami undergoes a vertical elongation. 
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FIG. 5. (Color online) Main: Phase diagram of the origami 
configuration before instability as a function of membrane 
length, L, and thickness, h: no origami (a), usual origami 
(■), origami with dewetting only (o), with ridge curvature 
only (•), with ridge curvature and dewetting (x). Top inset: 
ferrofluid dewetting (see left arrow) and ridge curvature (see 
right arrow) before instability. Bottom inset: "Anomalous" 
final state after the instability when adhesion is involved. 



one has Di — Z?i < that leads to a decrease of the mag- 
netie energy, i.e. AE^ < 0. The variation of the gravita- 
tional energy during the instability is AEg ~ pYgAzmc 
with Azmc the variation of the height of the mass center. 
At the onset of the instability, lS.Em compensates AEg , 
i.e. AE„ 



AEg = 0, and therefore 
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Az,nc = — / AM{H)dH, 
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For different sizes, L, of the membrane, Az,„c is de- 
duced from the measurements of z,„c before and after 
the instability occurring at Be- The origami aspect ratio 
a/b is also measured in the same way as for the ferrofluid 
drop alone on a plane (see above). For such a semi- 
ellipsoid, no theoretical expression for the demagnetiza- 
tion coefficient D is known. However, substituting the 
measured value a/b of the scmi-ellipsoid-like origami into 
the above expression of D{a/b) for a full ellipsoid leads 
to a measured value of the demagnetization coefficient 
Dm- The quantity Jg^" AMm{H)dHc is experimentally 
found as follows (m denoting a measured quantity). For 
each He = Sc/mo corresponding to each value of L, one 
computes AMm{H) by using the experimental values of 
Dm before and after the instability, by numerically solv- 
ing the above implicit equations given Ha and iJt,, and 
finally by iterating Hq from to He- Fig. |4]then shows 
Azmc as a function of J^^" AMm{H)dH(. for different val- 
ues of L. Both quantities are found proportional, and 
with a constant fi()/{2pg), i.e. half the coefficient pre- 
dicted in Eq. ([T]). Since magnetic energy mostly focuses 



near drop regions of small curvature (tip-like effect), the 
value of the magnetic energy for a semi-ellipsoid can be 
assumed to be half of the one of the full ellipsoid, i.e. 
AM = AMm/2. The experimental curve is therefore in 
good agreement with the predictions of Eq. ^ with no 
adjustable parameter. Finally, Azmc is measured to be 
independent of the thickness h at fixed L, emphasizing 
that elasticity does not play a significant role during this 
overturning instability. 

Capillarity and elasticity are at the origin of the capil- 
lary origami formation but are not involved in the thresh- 
old of the gravity-magnetic instability reported here. 
However, capillarity and elasticity effects can appear in 
the origami configuration before the instability: dewet- 
ting of the drop in the vicinity of the strongest curvature 
of the membrane (see left arrow in the top inset of Fig. 
[5]), as well as curvature of a ridge of the pyramid can 
occur (see right arrow). Both phenomena result from 
elasticity and capillarity since they depend on the mem- 
brane thickness h and length L according to the phase 
diagram shown in Fig. [5] For thin enough membranes 
and/or long enough ones, the elasticity energy ^ h^L^ 
is much smaller than the capillary ~ hPL^ and magnetic 
~ h^L^ ones. It thus leads to very small curvature radii 
of the membrane that rigidify the origami, and conse- 
quently favours the ferrofiuid dewetting (see zone 3 in 
Fig. [5)). For thinner h and/or longer L, as elasticity 
becomes negligible with respect to magnetic effects, the 
membrane bends inward and follows the elongated fer- 
rofluid shape, leading to the curvature of the pyramid 
ridges (see zone 4). These additional effects of elastic- 
ity and capillarity on the overturning instability deserves 
further studies. Note also that if an adhesion force exists 
between the substrate and the membrane, the instabil- 
ity leads to an anomalous final state: one corner of the 
membrane points up while the two others wrap the drop 
at its base (see bottom inset of Fig. [S]). 

Finally, the instability reported here could have some 
interest for potential applications (see Ref. [H). Indeed, 
the overturning instability is controlled in a non intrusive 
way and its scaling law is suitable for miniaturization as 
lower critical magnetic fields are required for smaller size 
membranes. 
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